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ABSTRACT 

We have used the Hubble Space Telescope Wide-Field and Planetary Camera 2, in combination 
with ground-based spectroscopy, to measure the integrated flux of galaxies at optical wavelengths — 
the extragalactic background light (EBL). We have also computed the integrated light from individual 
galaxy counts in the images used to measure the EBL and in the Hubble Deep Field. We find that flux 
in galaxies as measured by standard galaxy photometry methods has generally been underestimated by 
about 50%, resulting from missed flux in the outer, lower surface-brightness parts of galaxies and from 
associated errors in the estimated sky level. Comparing the corrected, integrated flux from individual 
galaxies with our total EBL measurement, we hnd that there is yet further light that contributes to the 
background that is not represented by galaxy counts, and that the total flux in individually detected 
sources is a factor of 2 to 3 less than the EBL from 8000 to 3000 A. We show that a significant fraction 
of the EBL may come from normal galaxies at z < 4, which are simply undetectable as a result of 
K-corrections and cosmological surface brightness dimming. This result is consistent with results from 
recent redshift surveys at z < 4. In the context of some simple models, we discuss the constraints placed 
by the EBL on evolution in the luminosity density at z > 1; while significant flux comes from galaxies 
beyond the current detection limits, this evolution cannot be tightly constrained by our data. 

Based on our measurements of the optical EBL, combined with previously published measurements of 
the UV and IR EBL, we estimate that the total EBL from 0.1-1000 /rm is 100±20 nW m“^ sr“^. If the 
total EBL were produced entirely by stellar nucleosynthesis, then we estimate that the total baryonic 
mass processed through stars is = 0.0062(±0.0022)/i“^ in units of the critical density. For currently 
favored values of the baryon density, Hb, this corresponds to 0.33 ± 0.12f2B. This estimate is smaller 
by roughly 7% if we allow for a contribution of 7/io.7 nW m“^ sr“^ to the total EBL from accretion 
onto central black holes. This estimate of suggests that the universe has been enriched to a total 
metal mass of O.21(±O.13)Z0 Db, which is consistent with other observational estimates of the cumulative 
metal mass fraction of stars, stellar remnants, and the intracluster medium of galaxy clusters in the local 
universe. 

Subject headings: Diffuse radiation — cosmology: observations — galaxies: evolution — galaxies: 
photometry 


1. INTRODUCTION 

The integrated optical flux from all extragalactic sources 
is a record of the stellar nucleosynthesis in the universe 
and the chemical evolution which has resulted from it. 
In Bernstein, Freedman, & Madore 2001a (henceforth. 
Paper I), we presented detections of the optical EBL in 
the HST/WFPC2 wide-band filters F300W (U300 , Aq ^ 
3OOOA), F555W (U555, Ao ~ 55OOA), and F814W (/gw, 
Aq ^ 8OOOA) based on simultaneous data sets from Hubble 
Space Telescope (HST) and Las Campanas Observatory 
(LCO). In Bernstein, Freedman, & Madore 200Ib (hence¬ 
forth, Paper H), we presented details of a measurement of 
the diffuse foreground zodiacal light which we use in Pa¬ 
per 1. Here we briefly summarize the results of Papers I 
and H and discuss the cosmological implications of these 
detections of the EBL. 

The majority of the EBL at UV to IR wavelengths is pro¬ 
duced by stars at restframe wavelengths of O.l-lO^m. Due 


to cosmic expansion, the EBL at U300, V555, and /814 po¬ 
tentially includes redshifted light from stellar populations 
out to z ~ 8 (the redshifted Lyman-limit cut-off of the 
/814 filter). Although stars themselves do not emit much 
light at wavelengths longer than lO^m, a complete census 
of the energy produced by stellar nucleosynthesis in the 
universe must consider the EBL over the full wavelength 
range 0.1-1000/im because dust in the emitting galaxies 
will absorb and re-radiate starlight, redistributing energy 
from nucleosynthesis into the thermal IR. 

With 8m-class telescopes and HST, the limits of resolved- 
source methods (i.e., number counts, redshift surveys, 
QSO absorption lines, etc.) for studying star formation 
in the universe are being extended to ever fainter levels; 
however, a direct measurement of the EBL remains an in¬ 
valuable complement to these methods. Galaxies with low 
apparent surface brightness — both intrinsically low sur¬ 
face brightness galaxies at low redshift and normal surface 
brightness galaxies at high redshift — are easily missed in 
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surface-brightness-limited galaxy counts and consequen- 
tally in follow-up redshift surveys. Identihcation, not to 
mention photometry, of faint galaxies becomes very uncer¬ 
tain near the detection limits. Even efforts to understand 
galaxy evolution, chemical enrichment, and star formation 
through QSO absorption line studies appear to be biased 
against chemically enriched, dustier systems, as these sys¬ 
tems can obscure QSOs which might lie behind them (Fall 
& Pei 1989, Pei & Fall 1995, Pettini etal. 1999). In con¬ 
trast, a direct measurement of the spectral energy distri¬ 
bution (SED) of the EBL from the UV to the far-IR is a 
complete record of the energy produced by star formation 
and is immune to surface brightness selection effects. 

In addition to energy originating from stellar nucleosyn¬ 
thesis, the EBL includes energy emitted by accreting black 
holes in quasars and active galactic nuclei. However, at op¬ 
tical wavelengths, the quasar luminosity functions at red- 
shifts z ^ 5 indicate that the optical luminosity density 
from quasars is a small fraction (~ 2.5%) of the that from 
galaxies (e.g. Boyle & Terlevich 1998). In addition, our 
measurement of the EBL excludes any point-like sources 
(of which there are 3 in our images), under the prior as¬ 
sumption that those sources are Galactic foreground stars. 
We therefore expect quasars to be a negligible source of 
flux in our measurements of the optical EBL. 

The contribution from active galactic nuclei (AGN) is 
more difficult to assess, as recent dynamical evidence 
(Richstone etal. 1998) indicates that massive black holes 
reside in most galaxies and sensitive optical spectroscopy 
(Ho etal. 1997a, 1997b) indicates that AGN have at least 
a weak contribution to more than 50% of nearby galax¬ 
ies. Nonetheless, simple accretion models, the total X-ray 
background, and the X-ray to far-IR spectral energy dis¬ 
tribution of AGN and quasars all indicate that the total 
contribution to the bolometric EBL fro m accretion onto 
central black holes is < 15% (see §^), and is emitted 
at thermal IR wavelengths. In principle, measurements 
of the EBL also constrain possible the total energy out¬ 
put from more exotic sources, such as gravitationally col¬ 
lapsing systems, brown dwarfs, and decaying particles (see 
Carr, Bond, & Hogan 1986, 1991 and Dwek etal. 1998 for 
discussions). 

The outline of the paper is as follows. In §^, we give an 
overview of the observations and methods used to measure 
the EBL as discussed in Papers I and H. In we sum¬ 
marize the individual measurements and associated errors 
we have obtained from each data set and the final EBL 
detections which result from them. In we compare 
the measured EBL with the integrated optical flux from 
resolvable sources as quantified by number counts and lu¬ 
minosity functions. In we quantify the contributions 
to the optical EBL which one might expected from sources 
which fall below the detection limits of the HDF based on 
explicit assumptions regarding the surface brightness, lu¬ 
minosity, and redshift distribution of galaxy populations 
in the universe. In §|[ we discuss models of the SED of the 
EBL based on these and recent results in the far infrared. 
Finally, in §0we discuss the total star formation and chem¬ 
ical enrichment history of the universe required to produce 
the bolometric flux of the EBL, and compare the inferred 
values to other observations of the total baryon fraction 
in stars and the metal mass density in the local universe. 


We abbreviate the adopted units ergs s ^ cm ^ sr ^ A ^ 
as cgs throughout. 

2. SUMMARY OF OBSERVATIONS 

As is true of all background measurements, the difficulty 
in measuring the optical EBL is in differentiating it from 
the much brighter foregrounds: terrestrial airglow, zodia¬ 
cal light (ZL), and diffuse Galactic light (DGL). Relative 
to the EBL flux at ~5000A, airglow and ZL are each more 
than 100 times brighter than the EBL. Along the most 
favorable lines of sight, the DGL is roughly equal in flux 
to the EBL. We have measured the EBL in a held which 
is out of the ecliptic plane and near the Galactic pole in 
order to optimally minimize the contributions of zodiacal 
light, DGL, and nearby stars (see Paper I). 

In the EBL measurement presented in Paper I, we 
have used three simultaneous data sets to isolate the dif¬ 
fuse EBL from the foreground sources: (I) absolute sur¬ 
face photometry taken with WFPG2 aboard HST using 
the wide-band Liters F300W (G300), F555W (V555), and 
F814W (/ 814 ); (2) low resolution (^ 300A) surface spec¬ 
trophotometry at 4000— 7 OOOA taken with the EOS, also 
aboard HST; and (3) moderate resolution (~2A) sur¬ 
face spectrophotometry taken with the Boiler and Chivens 
spectrograph on the 2.5m duPont telescope at LCO. We 
use the two HST data sets to measure the total mean hux 
of the night sky, including ZL, DGL and the EBL. We 
avoid terrestrial airglow all together by using HST for this 
measurement. We then use the LCO spectra to measure 
the absolute surface brightness of ZL in the same held and 
on the same nights as the HST observations. Finally, we 
estimate the small DGL contribution using a scattering 
model which is in good agreement with the observations. 
We then subtract the measured ZL and the modeled DGL 
from the total hux measured with HST/WFPC2 through 
each hlter and with HST/FOS. These measurements are 
described in detail in Papers I and H. Below, we summa¬ 
rize the observations, results, and accuracy of the individ¬ 
ual measurements which contribute to the EBL detection 
(see Table 

Bright galaxies brighter are not statistically well sam¬ 
pled in the 4.4 arcmin^ WFPC2 held of view. We have, 
therefore, masked out any sources brighter than V555 = 
23AB mag in the WFPC2 images before we measured the 
total sky hux. To do so, we used masks which are de¬ 
rived from the F555W images and extend to four times 
the isophotal radius in those data. We use the abbrevia¬ 
tion EBL23 as a reminder of this bright magnitude cut¬ 
off. The EBL23 detections can be combined with ground- 
based counts at V555 < 23AB mag to obtain the total EBL. 
The WFPC2 surface brightness measurements have ran¬ 
dom errors of < 1% and systematic uncertainties of I- 
2% of the total background hux. From the HST/WFPC2 
data alone, we can also identify a minimum hux from de¬ 
tectable sources. This minimum is given in Table |^, and 
the method used to obtain this result is summarized in 

The EOS spectra also provide a measurement of total 
hux. The random error per resolution element is around 
2.1%, and the systematic uncertainty over the full range 
is 3.5%. The ^ 14 arcsec^ EOS held of view and ^ 4% 
systematic uncertainty make the EOS less useful than the 
WFPC2 for measuring the EBL. However, most of sys- 
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tematic uncertainty is due to the poorly constrained solid 
angle of the aperture and aperture correction. Both of 
these are wavelength-independent errors, so that the FOS 
spectra do provide a ± 1 % measurement of the color of the 
total background, which is dominated by zodiacal light. 

The scattering which produces the ZL is well described 
by classical Mie theory for the large (> lO^m), rough dust 
grains which populate the zodiacal dust cloud. The scat¬ 
tering efficiency of the dust is only weakly wavelength de¬ 
pendent, so that the solar spectral features are well pre¬ 
served in the scattered ZL spectrum. However the broad 
band spectrum of the zodiacal light is redder than the solar 
spectrum by about 5% per lOOOA (see Paper II for further 
discussion) due to surface roughness of the grains, which 
decreases scattering efficiency at shorter wavelengths. The 
mean ZL flux in a narrow band can thus be measured from 
the apparent equivalent width of the solar Fraunhofer lines 
evident in its spectrum. Small color corrections can then 
be used to infer the full spectrum relative to that measure¬ 
ment. This requires moderate resolution spectra (~ 2A) 
with excellent flux calibration (± 1 %), which can only be 
obtained with ground-based observations, and then only 
at wavelengths where atmospheric emission lines are rela¬ 
tively weak. We have, therefore, measured the ZL in the 
range 4000-5100A using spectra taken at LCO. The re¬ 
sulting measurement has a statistical error of < 1 % and a 
systematic uncertainty of ^ 1.2%. This measurement has 
been extrapolated the 3000A and SOOOA WFPC2 band- 
passes using measurements of the color of the ZL from the 
FOS and ground-based LCO data. 

Within the Galaxy, there is both resolved flux from dis¬ 
crete stars and diffuse light (diffuse Galactic light, DGL) 
from starlight scattered by interstellar dust. Discrete stars 
can simply be resolved and subtracted in the WFPC2 im¬ 
ages. The intensity of the dust-scattered optical DGL and 
the 100 /rm thermal emission from the same dust are both 
proportional to the dust column density and the strength 
of the interstellar radiation field. To minimize the opti¬ 
cal DGL, our field was selected to have very low 100p,m 
emission. The remaining low-level DGL which does con¬ 
tribute has been estimated using a simple scattering model 
based on the dust column density and interstellar radi¬ 
ation field along the line of sight and empirical scatter¬ 
ing characteristics for interstellar dust. The predictions 
of this model are in good agreement with observations of 
the DGL from 2500-9000A (see Witt et al. 1997 and ref¬ 
erences therein). Finally, although isotropic line emission 
from warm interstellar gas is measured at all Galactic lat¬ 
itudes, the strongest line, Ha, does not lie within any our 
HST/WFPG2 bandpasses. The next strongest lines, [Oiii], 
are twenty times weaker and contribute negligibly to our 
results. 

The EBL cannot be measured in typical HST data. Our 
HST observations were scheduled to avoid contaminating 
scattered light from all anticipated sources: the bright 
Earth limb, the Moon, and off-axis stars. Also, obser¬ 
vations from LGO and HST were strictly simultaneous to 
guarantee that the ZL measured from the ground is ex¬ 
actly the contribution seen by HST. As an additional safe¬ 
guard, observations were scheduled in 3 visits, allowing us 
to look for possible off-axis scattered light with the satel¬ 
lite oriented at different roll angles, to safeguard against 
unidentified photometric anomalies with the instruments. 



4000 6000 8000 

Wavelength(A) 


Fig. 1. — Summary of EBL23 measurements, repeated from Paper 
I. Filled circles show the EBL23 obtained from surface photometry 
of the total background measured from HST/WFPC2, the zodia¬ 
cal light as measured from Las Campanas Obsermtory, and models 
of the diffuse galactic light as summarized in §H. The solid, dot¬ 
ted, and long-dashed error bars show the combined, systematic, and 
random la errors, respectively. The hatch—marked region shows the 
la uncertainty in the detected EBL due to uncertainty in ZL color. 
The lower limit arrows connected by a dashed line indicate the to¬ 
tal flux from individually photometered galaxies with magnitudes 
23 < F 555 < 30 AB mag in the HDF catalog. The u-shaped lower 
limit arrows show minEBL23, which is the flux as determined by 
ensemble photometry from galaxies with 23 < V 555 < 28 AB mag in 
the EBL fields. 


and to confirm the expected modulation in the ZL with 
the Earth’s orbital position. 

3. SUMMARY OF EBL DETECTIONS 

The individual measurements which are combined to ob¬ 
tain our detections of the EBL are summarized in Table 
^ We summarize our confidence intervals on the detected 
EBL23 in Table || and Figure |^. For comparison with 
the EBL23 fluxes, we have included in Table ^ the in¬ 
tegrated flux from individually photometered sources in 
the HDF, as measured using the photometry package FO- 
GAS (Jarvis & Tyson 1981, Valdes 1982) and published in 
Williams etai (1996). As the values in this table show, 
the mean EBL23 detections in each bandpass are more 
than 5x higher than the integrated flux in HDF galaxies 
as measured by standard photometry. 

To help understand this large difference between the 
detected EBL and the flux in HDF number counts, we 
have measured the total flux from resolved galaxies in our 
WFPG2 images (23 < V 555 < 27.5) using a method which 
we call “ensemble aperture photometry.” This method is 
uniquely suited to both our goal of a very accurate mea¬ 
surement of the ensemble flux of all galaxies in our images 
and to our data set, which has zero-point calibration ac¬ 
curate to ± 0 . 1 % over each image and negligible scattered 
light. This method is described in detail in Paper I and 
summarized below. 

Briefly, we identify the total flux from detectable galax- 
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Table 1 

Summary of Measurements 


Source 

Bandpass 

Data source 

Flux 

Random 

Systematic 

Total 

F300W 

hst/wfpc2 

33.5 

(± 4 . 9 %) 

[± 5.6%] 

Background 

F555W 

hst/wfpc2 

105.7 

(± 0.3%) 

[± 1 . 4 %] 


F814W 

hst/wfpc2 

72.4 

(± 0.2%) 

[± 1 . 4 %] 


F555W“ 

hst/fos 

111.5 

(± 0.7%) 

[± 2.8%] 

Zodiacal 

4600-4700A 

LCO 

109.4 

(± 0.6%) 

[± i.i%] 

Light 

F300W 

LCo'^ 

28.5 

(± 0.6%) 

[- 1 . 1 %,± 1 . 2 %] 


F555W 

LCO*^ 

102.2 

(± 0.6%) 

[- 1 . 1 %,± 1 . 1 %] 


F814W 

LCO*^ 

69.4 

(± 0.6%) 

[- 1 . 3 %,± 1 . 1 %] 

Diffuse 

F300W 

DGL model 

1.0 


[±25%,-50%] 

Galactic 

F555W 

DGL model 

0.8 


[±25%,-50%] 

Light 

F814W 

DGL model 

0.8 


[±25%,-50%] 


Note. —All fluxes are in units of 1x10 ®ergs s ^ cm ^ sr ^ A For a source with constant 
flux in F\, filters F300W, F555W, and F814W have effective wavelengths Ao(AA) =3000(700), 
5500(1200), and 8000(1500)A. Individual sources of error contributing to each measurement are 
summarized in Tables 3 & 4 of Paper I and Table 1 of Paper II. 

“Observed FOS spectrum, convolved with the WFPC2/F555W bandpass to allow direct com¬ 
parison with the WFPC2 results. 

^LCO measurement of zodiacal light have been extrapolated to the WFPC2 bandpass by apply¬ 
ing a correction for changing zodiacal light color with wavelength relative to the solar spectrum. 
The zodiacal light flux through the WFPC2 bandpasses was identified using SYNPHOT models, 
the uncertainty due to which is included in the uncertainty for the Alter calibration and is shared 
with the systematic uncertainty for the total background flux. 


Table 2 

EBL Results and Uncertainties 


Bandpass 

Random 

Systematic 

Combined 

EBL(±lo-) 

an (68%) 

as (68%) 

a (68%) 


Detected EBL23 (WFPC2 ± LCO)“ 




F300W 

2.1 

1.5 

2.5 

4.0 (±2.5) 

F555W 

0.6 

1.3 

1.4 

2.7 (±1.4) 

F814W 

0.4 

0.9 

0.0 

2.2 (±1.0) 

Minimum EBL (WFPC2) 

a 



F300W 

0.19 

0.13 

0.22 

3.2 (±0.22) 

F555W 

0.003 

0.009 

0.01 

0.89 (±0.01) 

F814W 

0.002 

0.007 

0.01 

0.76 (±0.01) 

Detected EBL23 (FOS + 

LGO)“ 



F555W 

0.7 

2.7 

2.8 

8.5 (±5.6) 


Flux from detected sources in HDF (m > 23 AB mag) 


F300W 

F450W 

F606W 

F814W 


0.66 

0.51 

0.40 

0.27 


Published number counts ^ 

F300W (18 < f/300 < 23 AB mag) 
F555W (15 < V555 < 23 AB mag) 
F814W (13 < /814 < 23 AB mag) 


0.27 (±0.05) 
0.49 (±0.10) 
0.65 (±0.13) 


Note. —All fluxes and errors are given in units of 10 ® ergs s ^ cm ^ sr ^ A 

“The systematic and statistical errors have been combined assuming a flat or Gaussian probabil¬ 
ity distribution, respectively, as discussed in Paper I. We equate Icr combined errors with the 68% 
confidence interval, as the combined errors are nearly Gaussian distributed. Individual sources of 
error contributing to these totals are summarized in Tables 3 and 4 of Paper I and in Table 1 of 
Paper II. 

^Estimated errors correspond to uncertainties in the fits to published galaxy counts. The values 
given correspond to 0.081 X 10“^®, 0.46 X 10“^'^, and 1.5 X 10“^'^ in units of ergs s“^ cm“^ sr“^ 
Hz“^ and are consistent with those used in Pozzetti et al. (1998). 
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ies fainter than i4ut = 23 AB mag by simply masking out 
galaxies with V < ABmag (and all stars) and av¬ 
eraging the flux of every pixel in the frame. From this, 
we obtain the mean surface brightness of foregrounds plus 
all extragalactic sources, or the average surface brightness 
per pixel from “objects -I- sky.” We then mask out all 
detected objects using standard detection apertures (twice 
the isophotal area) and calculate the average flux in the 
remaining pixels. From this, we obtain the mean surface 
brightness outside of the galaxy masks, or the average sur¬ 
face brightness per pixel from “sky.” The difference be¬ 
tween these two measurements is then the ensemble sur¬ 
face brightness of all sources within the area of the masks. 
This assumes that the sky level is uniform, which is the 
case to better than 1% accuracy in our images. By varying 
the bright magnitude cut-off (14ut) we choose for measur¬ 
ing “objects -I- sky,” we can isolate the flux coming from 
sources fainter than t 4 ut- 

As discussed in Paper I, we find that the recovered flux 
increases steadily with increasing mask size. For exam¬ 
ple, roughly 20% of the light from galaxies 4 magnitudes 
above the detection limit lies at radii •\/2?’iso < r < 4riso 
(see Figure H), beyond the standard-size galaxy apertures 


(v 2 riso) used in faint galaxy photometry packages, such as 
SExtractor (Bertin & Arnouts 1996) or FOCAS. Because 
estimates of the sky level in standard photometry packages 
come from just beyond the detection apertures, these sky 
estimates will include some fraction of the galaxies’ light 
and will doubly compound this error. In addition, because 
galaxy apertures start to significantly overlap in our im¬ 
ages and the HDF images when they extend to r ^ 4riso, 
we find that some flux from the wings of detected galaxies 
will inevitably contribute a pedestal level to the mean sky 
level in any image. We have estimated this pedestal level 
by Monte Carlo simulations as described in Appendix B 
of Paper I. The pedestal is independent of field, but does 
dependent on the detection limits and surface brightness 
characteristics of the data. For Veoe HDF images, this 
pedestal level is roughly 10% of the total flux from V > 23 
AB mag galaxies and, again, this error is compounded by 
the fact that any flux at radii beyond galaxy apertures 
can be include in the “sky” estimate. The true flux from 
V > 23 AB mag galaxies in the HDF is therefore alm ost 
twice what is recovered by standard methods (see §4.1). 

Using different values of Ucut, we can quantify system¬ 
atic errors in faint galaxy photometry as a function of the 
isophotal surface brightness limit of the data, fj-iso, and the 
central surface brightness of the source, /ig. The smaller 
the value of = /Xiso — is for a particular galaxy, the 
larger the photometric error in standard aperture photom¬ 
etry. This problem has been discussed in the literature at 
length with respect to low surface brightness galaxies at 
low redshifts based on extrapolations of simple exponen¬ 
tial light profiles (Disney 1976, Disney & Phillips 1983, 
Davies 1990); the same principle begins to apply to nor¬ 
mal surface brightness galaxies which have low apparent 
surface brightness at higher redshifts due to (1 -I- sur¬ 
face brightness dimming (Dalcanton 1998). 

Finally, we note that the direct measurements of the 
EBL23 in our data — based on surface photometry of 
the total integrated background, zodiacal light, and dif¬ 
fuse galactic light — are 2cr — 3a detections. However, 


the mean flux from detected sources is obviously an abso¬ 
lute minimum value for the EBL. Therefore, the strongest 
lower limit we can place on the flux from sources fainter 
than U = 23 AB mag (EBL23) is the mean flux in detected 
galaxies as measured by ensemble aperture photometry in 
our WFPC2 data and the HDF. The strongest upper lim¬ 
its we can place on EBL23 are the 2a upper limits of our 
direct measurements of the EBL23. In Table ^ we list (1) 
our direct measurements of the EBL23, and (2) the min¬ 
imum values for the EBL23 (minEBL23) as identified by 
ensemble aperture photometry of detected sources in our 
WFPC2 data and the HDF. For comparison, the flux in 
published HDF galaxy counts and ground based counts 
are also listed there. 

4. COMPARISON WITH NUMBER COUNTS AND 
LUMINOSITY FUNCTIONS 

Whether the light originates from stellar nucleosynthe¬ 
sis, accretion onto compact objects, or gravitationally col¬ 
lapsing stellar systems, the total optical flux escaping from 
detected galaxies is quantified by number counts and lumi¬ 
nosity functions. To the detection limits, number counts 
and luminosity functions contain exactly the same infor¬ 
mation regarding the integrated background light: the in¬ 
tegrated flux from resolved sources is the same whether or 
not you know the redshift of the sources. However, in the 
context of predicting the EBL flux, luminosity functions 
contain information about the intrinsic flux distribution 
of the sources and thus allow us to estimate the flux from 
sources beyond the detection limits with better defined 
assumptions. In the following sections, we compare our 
EBL detections with the integrated flux obtained by both 
methods. Dust obscuration in the emitting sources will 
clearly reduce the UV and optical flux which escapes, but 
the EBL, number counts, and luminosity functions are all 
measurements of the escaping flux; the relative compar¬ 
isons discussed in this section are therefore independent of 
dust extinction. 


4.1. Number Counts 

Using “ensemble aperture photometry” to measure the 
total flux from galaxies as a function of magnitude in our 
V 555 and Jgi 4 images of the EBL field, we find that the 
standard photometry methods used to produce the HDF 
catalog systematically underestimate the flux from each 
source, as summarized in (see §10 and Appendix B of 
Paper I for a thorough discussion). We use these results 
to derive flux corrections as a function of ApL = piiso — /ig 
(isophotal minus central surface brightness), which are es¬ 
sentially aperture corrections. These aperture corrections 
are similar to those found by other authors (c.f. Small et al. 
1995) and are a natural consequence of integrating an ex¬ 
tended light profile to an insufficient radius. This effect 
can be quantified for exponential or de Vaucouleurs pro¬ 
files, as in Dalcanton (1998). However, the corrections we 
show here are empirical measurements and assume nothing 
about the light profiles of the sources. 

The corrections we derive for the two bandpasses (see 
Figure]^) are very similar functions of A^, which indicates 
that the profiles of detected galaxies are not a strong func¬ 
tion of wavelength over the baseline of observed V to I. 
However, we note that a particular value of A/i occurs 
at a brighter AB magnitude in Jgi 4 than in Vggg because 
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Fig. 2.— Aperture corrections as a function of A/i = fiiso — fio 
(isophotal minus central surface brightness) derived by “ensemble 
aperture photometry” of the EBL field for V (thick line) and I (thin 
line). The mean apparent magnitude in Veoe Isi 4 corresponding 
to a particular value of Afi in the HDF images is indicated by the 
x-axes at the top of the plot. Error bars show the Icr statistical error 
in the mean corrections derived from 18 WFPC2 images of the EBL 
field (six exposures, three WF chips). 


the limiting isophotal level (sky noise) in J 814 is 0.6 AB 
mag brighter than in Veoe- The corrections are therefore 
larger in /814 than they are at the same AB magnitude 
at keoe- The corrections in both bands include a correc¬ 
tion which compensates for overestimates in the sky flux 
from foreground sources (the pedestal sky level described 
in §^) . This correction, which accounts for errors in the lo¬ 
cal sky estimate, ranges from 0.1-0.3 mag, monatonically 
increasing towards fainter magnitudes. As in Vqqq and 
/ 814 , aperture corrections for U 300 band will depend on 
the prohles of galaxies at C /300 and the surface brightness 
limits of the data. However, the very low signal-to-noise 
ratio of our F300W images prevents us from determining 
aperture corrections in that bandpass. The C /300 photom¬ 
etry is discussed further below. 

We have applied the aperture corrections we derive to 
the individual objects in the HDF Vooo and /814 catalogs 
(Williams et al. 1996), which fractionally increases the flux 
of each galaxy. For example, while galaxies in the HDF 
catalog with Vooo 30 AB mag have well-detected cores, 
less than 30% of their total flux is recovered: the total 
flux of a galaxy measured to have Vooo 30 ABmag by 
standard photometry methods is actually closer to Vooo ~ 
28 ABmag. The corrected and uncorrected (raw) galaxy 
counts and corresponding integrated flux with magnitude 
are compared in Figures ^ and The integrated flux 
of the corrected galaxy counts roughly corresponds to the 
minimum value of EBL23, as the aperture corrections were 
derived from the calculation of the minimum EBL23 in our 
data. Statistical variations in galaxy counts between helds 
are to be expected. 

The aperture corrections we apply clearly have a sig¬ 
nificant impact on the slope of faint number counts. To 



Veoa AB mag 



VgQg AB mag 


Fig. 3.— The upper panel shows the galaxy counts as published in 
the HDF catalog (filled circles) with \/N error b ars a nd the corrected 
number counts (open circles), as described in §4.1. The solid lines 
show fits to the raw number counts, which change slope slightly 
around Vqoq = 26. The fit to the corrected counts is indicated by 
the dashed line to the detection limit and a dotted line beyond. No 
change in slope is apparent at the faint end for the corrected counts. 
Ail slopes are given in the text. In the lower panel, we plot the 
integrated flux corresponding to the galaxy counts with the same line 
types as in the upper panel. The data point and Icr error bar mark 
the value of EBL23 (converted to Veoe from the V 555 band). The 
corresponding ±lcr error range is emphasized by the hatch-marked 
region. For comparison, the lower limit arrow shows 2cr lower limit of 
minEBL23, the integrated flux from detected sources with Vgoe > 23 
AB mag. 




Ig ,4 AB mag 


Fig. 4.— The same as Figure ^ but for the 78i4—band. As for 
the Veoe^band, the raw /814 counts show a slight change in slope 
around 24—26 AB mag, while the corrected counts do not. All slopes 
are given in the text. The integrated flux of the raw and corrected 
counts are compared to our EBL23 7814 detections in the lower panel, 
as in Figure ^ 
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_ Fig. 5.— The raw and corrected number counts from Figure 
ra compared to number counts from the literature, labeled by first 
author. The lines indicate fits to the data using the relation N oc 
10“"*; a = 0.48 ±0.1 and a = 0.33 ± 0.01 at the bright and faint 
ends, respectively. Note that the slope of the corrected HDF counts 
is well matched to that at brighter magnitudes. 


quantify this, we fit both the raw and corrected number 
counts with the usual relationship between apparent mag¬ 
nitude and surface number density, N oc 10“™, where N is 
the number of galaxies per magnitude per square degree. 
For the raw Vgoe counts, we find that the data exhibit a 
change in slope around 24-26 AB mag. A single fit over 
the range 22 < Veoe < 29.5 AB mag gives a = 0.24 ± 0.01 
with a per degree of freedom (y^/dof) of 1.5. Fitting 
the counts brighter and fainter than 26 AB mag, respec¬ 
tively, we find ab = 0.28 ± 0.02 with y^/dof = 0.9 and 
a/ = 0.21 ± 0.01 with y^/dof = 1.2 (solid lines in the up¬ 
per panel of Figure §). We ascribe this change in slope to 
the onset of significant photometry errors. 

For the corrected Veoe, counts we find that the full 22 < 
Veoe < 27.5 AB mag range is well fit by a slope of a = 
0.33 ± 0.01 with y^/dof = 0.60 (dashed line in the upper 
panel of Figure y). This result suggests that photometry 
errors are responsible for the change in slope at the faint 
end of the HDF counts, and that N{m) does not, in fact, 
significantly decline before the detection limit of the HDF 
at Veoe- In addition, while the integrated flux in the raw 
galaxy counts has converged by the apparent detection 
limit of the HDF, the flux from the corrected galaxy counts 
has not (see the lower panel of Figure ||). 

We find similar results for the Isi 4 counts (see Fig¬ 
ure H). As for Veoe, the raw /814 counts display a slight 
change in slope around 24-26 AB mag. We find a slope of 
at = 0.25 ± 0.01 with y^/dof = 0.6, and a/ = 0.19 ± 0.02 
with y^/dof = 2.0, brighter and fainter than 26 AB mag, 
respectively. For the full range 22 < J 814 < 29.5 AB 
mag, we find a = 0.22 ± 0.01 with y^/dof = 2.1. For 
the corrected / 8 i 4 counts, we find a = 0.34 ±0.01 with 
y^/dof = 0.8 at 22 < J 814 < 27 AB mag. 

In Figure H, we show the raw and corrected HDF counts 
relative to 14- and i?-band counts available in the litera¬ 



bi 4 mag 

Fig. 6. — Same as Figure ^ but for the /—band counts. The lines 
indicate fits to the data using the relation N oc 10“™; a = 0.52±0.1 
and a = 0.34 ± 0.01 at the bright and faint ends, respectively. Note 
that slope of the corrected HDF counts is well matched to the slope 
at brighter magnitudes, and that the slope of the I and V band 
counts are similar at all magnitudes. 


ture for 14 > 15 AB mag. We have converted all of the 
published counts to 14-band AB mag by applying constant 
offsets consistent with those in Fukugita, Shimasaku, & 
Ichikawa (1995). These incorporate mean K-corrections 
based on the mean redshift corresponding to the apparent 
magnitude of the sample. Differences between filters will 
have some affect on the slope of counts in surveys which 
cover a large range of redshift (apparent magnitude) due 
to changing galaxy colors and K-corrections with increas¬ 
ing redshift, but these effects will average out between the 
multiple surveys shown. This plot shows that the aperture 
corrections we have applied to the HDF sources produce 
number counts which have a slope consistent with the slope 
found at brighter magnitudes. 

In Figure we show the same plot for the /-band. 
Again, the corrected / 8 i 4 counts display a slope which is 
similar to that found at magnitudes brighter than 23 AB 
mag. Note also that slope of the counts at < 25 AB mag in 
14 and / are the same to within the statistical errors. The 
aperture corrections we apply to the HDF counts at Vgoq 
and /814 extend this agreement to the current detection 
limits. The corrected counts imply that the faintest galax¬ 
ies detected do not exhibit a significantly steeper slope in 
VgQg than in / 8 i 4 , in contrast with the raw galaxy counts. 
This is an important constraint on galaxy evolution mod¬ 
els. 

Although the signal to noise in the t /300 data is too 
low to allow us to obtain accurate aperture corrections 
at that wavelength, the minimum EBL23 at C /300 implies 
consistent colors for faint and bright galaxies at C/ —14, as 
in 14 — / (see Figure and Table ^). We note, also, that 
the color of the integrated flux from galaxies is consistent 
with the color of the total background light within 2a. In 
other words, no exotic population of sources is required to 
produce the detected background. 




Detections of the Optical EBL (III): Cosmological Implications 


The lack of turnover in the corrected counts strongly 
suggests that sources do exist at apparent magnitudes be¬ 
yond the present detection limit. If we impose no limit 
on the apparent magnitude of sources and simply extrap¬ 
olate the galaxy counts beyond Veoe ~ 27.5 AB mag using 
a = 0.33 (dotted line in Figure ||), we obtain a prediction 
for the total integrated EBL23 of I.3xI0“® cgs, which is 
Icr below the measured value in the EBL field. In this case, 
the predicted EBL23 converges around Veoe 50 AB mag, 
which is signihcantly fainter than a dwarf galaxy at z ~ 6 . 
However, very little flux is obtained from the faintest bins. 
If we impose the limit Veoe ~ 38 AB mag as the faintest 
apparent magnitude for a realistic source (e.g., a dwarf 
galaxy with My ~ —10 AB mag at z ~ 4), we obtain a 
flux of I.2x I0“® cgs. The flux from sources with /814 > 23 
AB mag is 1.3x10“® cgs if we adopt a = 0.34, with the 
flux converging around /814 ~ 60 AB mag. Adopting a 
more realistic faint cut-off of ~ 38 AB mag, as discussed 
for Veoe, we obtain a total flux of 1 . 2 x 10 “® cgs, Icr below 
the mean detected value of EBL23 at /814 (see Figure ^). 

In order to obtain a cumulative flux equal to the mean 
detected EBL (or the upper limit) from sources brighter 
than ~ 38 mag, the slope of the galaxy counts in the 
range 28-38 AB mag would clearly need to increase at 
some point beyond the current detection limit. For exam¬ 
ple, the total flux from sources 23 < Veoe < 38 AB mag 
will produce the mean detected EBL if the sources with 
23 < Veoe < 28 AB mag obey a slope of a = 0.33 and 
sources with 28 < Veoe < 38 AB mag obey a = 0.42. We 
stress, however, that the total flux obtained from sources 
with such a steep faint-end slope is critically dependent on 
the cut of magnitude: the total flux reaches 5.1x10“® cgs 
if we integrate the counts to 50 AB mag, and to 9.Ox 10“® 
cgs if we integrate to 60 AB mag. Recall that our 2cr upper 
limit on EBL23 at Veoe is 5.0x10“® cgs. For a = 0.35 at 
E > 28 AB mag, the integrated flux reaches 1.37, 1.51, and 
1.57x10“® cgs (converged) for faint cut-off magnitudes of 
40, 60, and 80 AB mag respectively. Although it is ob¬ 
viously impossible to place firm constraints on the num¬ 
ber counts beyond the detection limit, as they may change 
slope at any magnitude, we conclude that it is very unlikely 
that the slope beyond Veoe ~ 28 AB mag is steeper than 
a = 0.42. If the slope continues at 0.33 < a < 0.35, the 
EBL23 reaches a roughly 1.3-1.5x10“® cgs by Eeoe ~ 40 
AB mag, < Icr below our detected value. 

Similarly, for the /-band the integrated flux from sources 
matches the mean detected EBL23 if the sources with 
23 < /814 < 27 AB mag obey a slope of a = 0.34 and 
sources with 27 < /814 < 39 AB mag obey a = 0.42. For 
those slopes, the total flux reaches the 2 a upper limit of the 
EBL23 at /814 by 50 AB mag. For a = 0.36, slightly above 
the slope we find for the corrected counts, the integrated 
flux reaches 1.31, 1.58, and 1.62x10“® cgs (converged) for 
faint cut-off magnitudes of 40, 60, and 80 AB mag, respec¬ 
tively. As for the E-band, we conclude that it is unlikely 
that the /-band faint-end slope is steeper than 0.42 at 
any magnitude. For a slope of 0.34 < a < 0.36 for / > 27 
AB mag, the EBL reaches 1.2-1.3x10“® by /814 ~ 40 AB 
mag, Icr below our detected value. 

In summary, we conclude from the corrected number 
counts shown in Figures l“l that sources are likely to 
exist beyond the detection limit of the HDF. Furthermore, 
if the number counts continue with the slope we measure 



Wavelength (A) 

Fig. 7. — The local luminosity density constructed by combin¬ 
ing the spectral energy distributions of E/SO, Sab, and Sc galaxies 
weighted according to the type=dependent luminosity functions as 
described in 34 .2| and Equation 


at the faintest levels, then the predicted EBL23 is within 
Icr of the detected EBL23 at both Eeoe and / 814 . If we 
extrapolate beyond the detection limits assuming the slope 
found from the corrected number counts, we find that less 
than 50% of EBL23 comes from sources beyond the current 
detection limit at Eeoe or /814 — the majority of the light 
contributing to EBL23 comes from sources which can be 
individually detected. 

Finally, we note that our ensemble photometry method 
yields a statistical correction for the light lost from the 
wings of galaxies beyond the detection isophote. This light 
cannot, by dehnition, be recovered by standard single¬ 
object photometry. In contrast, the ensemble photometry 
method effectively adds together the light beyond the de¬ 
tection isophote from many galaxies to enable a signihcant 
measurement of that light. 

4.2. Luminosity Functions 

In this section, we compare the detected EBL with 
the EBL predicted by luminosity functions measured as 
a function of redshift. To avoid unnecessary complications 
in defining apparent magnitude cut-offs, and to facilitate 
comparison with other models of the luminosity density 
as a function of redshift, we compare luminosity functions 
with the total EBL rather than with EBL23, as in the 
previous section. To do so, we combine the EBL23 flux 
measured in Paper I with the flux from number counts at 
brighter magnitudes, as given in Table ^ Sy stematic er¬ 
rors in photometry of the sort discussed in ^.l| are likely to 
be relatively small for redshift surveys because the objects 
selected for spectroscopic surveys are much brighter than 
the limits of the photometric catalogs (although see Dal- 
canton 1998 for discussion of the effects of small, system¬ 
atic photometry errors on inferred luminosity functions). 
We have not tried to compensate for such effects here. 

The integrated flux from galaxies at all redshifts is given 
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Fig. 8. — The upper panel shows the spectrum of the EBL calcu- 
^ted by integrating the luminosity density over redshift (Equation 
W) for constant luminosity density, passively evolving luminosity den¬ 
sity, and evolution of the form C(\^z) = £(A,0){1 -h with 

/Z{A, z) constant at 2: > 1. The lower panel shows S{X) for the three 
cases of (1 -has labeled in the figure and described in the text. 
Line types and hatch-marked regions in upper panel correspond to 
values of (5(A) with the same line type in the lower panel. Filled 
circles show the mean EBL detections with 2cr error bars. The error 
bars are dashed where they extend below the cumulative flux in de¬ 
tected sources — the minimum values for the EBL — indicated by 
the lower limit brackets. 


Fig. 9.— The three panels show the spectrum of the EBL calcu¬ 
lated assuming C{X, z) = C{X, 0)(l-h2:)^^'^^ over the range 0 < z < 1.5 
(bottom panel), 0 < z < 2 (middle), and 0 < 2 < 3 (top). In all 
cases, the luminosity density is held constant beyond the indicated 
redshift limit. The hatch-marked regions each show the dilcr range 
of CFRS values for 5(A), as in Figure fel—The integrated EBL as a 
function of redshift is shown in Figure^. Luminosity density as a 
function of redshift is shown in Figure h for some combinations of 
5(A) and the redshift cut-off for evolution. TJie filled circles, error 
bars and lower limit symbols are as in Figure R 


by 

/(A,0) = ^ (1) 

in which Vc{z) is the comoving volume element, Dl(z) is 
the luminosity distance, Aq is the observed wavelength, 
and Az = Ao(l + 2 :)“^ is the rest-frame wavelength at 
the redshift of emission. To compare the detected EBL 
to the observed luminosity density with redshift, £(A, z), 
we begin by constructing the SED of the local luminosity 
density as a linear combination of SEDs for E/SO, Sb, and 
Ir galaxies, weighted by their fractional contribution to the 
local B-band luminosity density: 

£(A,0) = ^A(B,0)^, (2) 

in which the subscript i denotes the galaxy Hubble type 
(E/SO, Sb, or Ir), /i(A) denotes the galaxy SED (the flux 
per unit rest-frame wavelength), and Ci{B,0) is the B- 
band, local luminosity density in ergs s“^ Mpc“^. To 
produce the integrated spectrum of the local galaxy pop¬ 
ulation, we use Hubble-type-dependent luminosity func¬ 
tions from Marzke etal. (1998) and SEDs for E, Sab, and 
Sc galaxies from Poggianti (1997). We adopt a local lu¬ 
minosity density of = 1.3 x 10 ®/iLqMpc“^, consistent 
with the Loveday etal. (1992) value adopted by CERS and 
also with Marzke etal. (1998)Jj The spectrum we obtain 
for £(A, 0) is shown in Eigure Q. 

^For a B-band solar irradiance of Lq = 4.8 X 10^®ergs 



Fig. 10.— The integrated EBL at V555 contributed as a function 
of increasing redshift from z = 0 to z = 10. As marked in the figure, 
the lines show the integrated flux for no evolution in the luminosity 
density, passive evolution, and evolution of the form C(X,z) oc (1 -h 
^)5C-^) for the — 1(T, mean, and H-lcr range of 5(A) values determined 
by CFRS. For each 5(A), we plot the growth of the EBL with redshift 
if C(X,z) is held constar^ at 2:_> 1, z > 1.5, z > 2, and 2 > 3, 
corresponding to Figures g and g. 


C{B,0) = 6.1 X 10^^h ergs s ^ Mpc ^ = 4.0 x lO^^/iso W Hz ^ 
Mpc“^. 
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Detections of the Optical EBL (III): Cosmological Implications 


We note that the recent measurement of the local lumi¬ 
nosity function by Blanton etal. (2001) indicates a factor 
of two higher local luminosity density than found by pre¬ 
vious authors. Previous results are generally consistent 
with Loveday et al. to within 40%. Blanton et al. attribute 
this increase to deeper photometry which recovers more 
flux from the low surface brightness wings of galaxies in 
the ir sa mple relative to previous surveys (see discussions 
in §4.1), and photometry which is unbiased as a function of 
redshift. For the no-evolution and passive evolution mod¬ 
els discussed below, the implications of the Blanton et al. 
results can be estimated by simply scaling the resulting 
EBL by the increase in the local luminosity density. Al¬ 
though the Blanton et al. results do not directly pertain 
to the luminosity functions measured by CFRS at red- 
shifts z > 0.2, they do suggest that redshift surveys at 
high redshifts will underestimate the luminosity density, 
as discussed by Dalcanton (1998). 

In the upper panel of Figure we compare the EBL flux 
we detect with EBL flux predicted by five different models 
for £(A, zL using the local luminosity density derived in 
Equation g as a starting point. For illustrative purposes, 
the first model we plot shows the EBL which results if we 
assume no evolution in the luminosity density with red¬ 
shift, i.e. £(A, z) = £(A, 0). The number counts them¬ 
selves rule out a non-evolving luminosity density, as has 
been discussed in the literature for over a decade; incon¬ 
sistency between the detected EBL and the no-evolution 
model is just as pronounced. The predicted EBL for the 
no-evolution model is a factor of 10 fainter than the de¬ 
tected values (filled circles). These are 1.7cr, 2.1cr, and 2.2cr 
discrepancies at C/300, 1455 and /814, respectively. More 
concretely, the no-evolution prediction is at least a fac¬ 
tor of 12X, 4x, and 3.7x lower than the flux in individ¬ 
ually resolved sources at C/300, 1455 and /814 (lower-limit 
arrows). Note that the no-evolution model still underpre¬ 
dicts the EBL if we rescale the local luminosity density to 
the Blanton etal. (2001) values. This model demonstrates 
the well-known fact that luminosity density is larger at 
higher redshifts. 

The second model we plot in Figure ^ shows the effect 
of passive evolution on the color of the predicted EBL. In 
this model, we have used the Poggianti (1997) SEDs for 
galaxies as a function of age for Hq = 50 km s“^ Mpc“^ 
and go = 0.225. In the Poggianti models, stellar popula¬ 
tions are 2.2 Gyrs old at a z ~ 3. The resulting £(A, z) is 
bluer than the no-evolution model due to a combination 
of K-corrections and increased UV flux for younger stellar 
populations. The passive evolution model does provide 
a better qualitatively match to the SED of the resolved 
sources (lower limits) and EBL detections (filled circles); 
however, it is still a factor of 3x less than the flux at 
C/ 300 , and a factor of 2x less than the flux we recover from 
resolved sources at 1455 and / 814 . For the adopted local lu¬ 
minosity density and Poggianti models, passive evolution 
is therefore not sufficient to produce the detected EBL. 
Again, the passive evolution adopted here still underpre¬ 
dicts the EBL if we rescale the local luminosity density 
by a factor of two to agree with the Blanton etal. (2001) 
value. 

As a fiducial model of evolving luminosity density, we 
adopt the form of evolution implied by the CFRS redshift 
survey (Lilly etal. 1996, hereafter CFRS) and Lyman-limit 


surveys of Steidel etal. (1999): C{\,z) = £(A, 0)(l-|-z)'^^'^^ 
over the range 0 < z < 1 and roughly constant luminos¬ 
ity density at 1 < z < 4. The remaining three models 
shown in Figure || test the strength of evolution of that 
form which is allowed by the EBL detections. The hatch- 
marked region shows the EBL predicted for values of (5(A) 
which represent the ± 1 ct range found by CFRS for the 
redshift range 0 < z < 1. The value of the exponent (5(A) 
is indicated in the lower panel of Figure ||, and the hatch- 
marked region reflects the uncertainty in the high redshift 
luminosity density due to the poorly constrained faint- 
end slope of the luminosity functions. This ± 1 (t range of 
the predicted EBL is consistent with the detected EBL at 
C/ 300 , but is inconsistent with the EBL detections at V555 
and /814 at the la level of both model and detections. It 
is, however, consistent with the integrated flux in detected 
sources at V555 and /814. 

To test the range of evolution allowed by the full ±2(T 
range of the EBL detections, we can explore two possibil¬ 
ities: (1) stronger evolution at 0 < z < 1, shown in Figure 
1 ^ and (^) evolution continuing beyond z = 1, shown in 
Figures ||, |l^, and Addressing the possibility of con¬ 
stant luminosity density at z > 1, the dashed line in the 
upper panel of Figure || shows the EBL predicted by the 2cr 
upper limit for (5(A) from CFRS; the dot-dashed line cor¬ 
responds to the value of (5(A) required to obtain the upper 
limits of the EBL detections at all wavelengths. Note that 
the latter implies a value for £(4400A, 1) which is ^ lOx 
higher than the value estimated by CFRS. This result em¬ 
phasizes that the 2a interval of the EBL detections span 
a factor of 4 in flux at 4400A, and thus the allowed range 
in the luminosity density for A < 4400A and 0 < z < 1 is 
similarly large. Also, for each model in which the luminos¬ 
ity density is constant at z > 1, less than 50% of the EBL 
will come from beyond z = 1 due to the combined effects 
of K-corrections and the decreasing volume element with 
increasing redshift (see Figure |T^ . 

Evolution continuing beyond z = 1 is possible if the 
Lyman-limit-selected surveys have not identified all of the 
star formation at high redshifts, and estimates of the lumi¬ 
nosity density at 3 < z < 4 are subsequently low. Figures 
1 ^ an(l 1 ^ show the EBL predicted by models in which the 
luminosity density increases as (1 -|- z)'^*^'*'^ to redshifts of 
z = 1.5, 2 and 3. Clearly, significant flux can come from 
z > 1 if the luminosity density continues to increase as 
a power law. The rest-frame C /300 luminosity density is 
plotted as a function of redshift in Figure |l^ for limiting 
values of the cut-off redshift for evolution and (5(A). Al¬ 
though the strongest evolution plotted over-predicts the 
detected EBL, our detections are clearly consistent with 
some of the intermediate values of the 5{X) and increas¬ 
ing luminosity density beyond z = 1. For example, the 
mean rate of increase in the luminosity density found by 
CFRS can continue to redshifts of roughly 2.5-3 without 
over-predicting the EBL. 

In all models, we have adopted the same cosmology 
{h = 0.5 and D = 1.0) as assumed by CFRS and Stei¬ 
del etal. (1999) in calculating £(A, z) and (5(A). Although 
the luminosity density inferred from these redshift surveys 
depends on the adopted cosmological model, the flux per 
redshift interval is a directly observed quantity. The EBL 
is therefore a directly observed quantity over the redshift 
range of the surveys, and is also model-independent. To 



Bernstein, Freedman, & Madore 


11 



Fig. 11 . — The luminosity density at f/300 fumution of redshift 
corresponding to limiting cases plotted in Figure ^ The hatch- 
marked region indicates the dilcr range given by the CFRS measure¬ 
ments of (5(A) over the range 0 < .z < 1 . The horizontal line segments 
show the luminosity density corresponding to the —la limit for (5(A) 
held constant at 2 > 1 . 5 ; the mean value for (5(A) held constant at 
z > 2\ and the -1-1(7 limit for (5(A) held constant at 2 > 3. 


the degree that the luminosity density becomes uncon¬ 
strained by observations at higher redshifts, the EBL does 
depend on the assumed (not measured) luminosity den¬ 
sity and on the adopted cosmology through the volume 
integral. Although dependence of the predicted EBL on 
Hq cancels out between the luminosity density, volume 
element, and distance in Equation Hq has some im¬ 
pact through cosmology-dependent time scales, which af¬ 
fect the evolution of stellar populations. If the luminosity 
density is assumed to be constant for z > 1 , the predicted 
EBL increases by 25% at V 555 for {Qm = 0.2, Ha = 0) and 
corresponding values of d(A), and decreases by 50% for 
(0.2,0. 8 ). The luminosity densities corresponding to the 
2cr upper limit of the detected EBL change by the same 
fractions for the different cosmologies if C is constant at 
z > 1. Similarly, for models in which the luminosity den¬ 
sity continues to grow at z > 1 , the luminosity density 
required to produce the EBL will be smaller if we adopt 
(0.2,0) than (1,0), and smaller still for (0.2,0.8). The exact 
ratios depend on rate of increase in the luminosity density. 

Several authors (Treyer etal. 1998, Cowie etal. 1999, 
and Sullivan etal. 2000) have found that the C at UV 
wavelengths (2000-2500A) is higher than claimed by CFRS 
(2800A) in the range 0 < z < 0.5 and have found weaker 
evolution in the UV luminosity density, corresponding to 
d(2000A) ~ 1.7. The implications for the predicted EBL 
can be estimated from the plots of the £{ 1 / 300 , z) shown 
in Figure and the corresponding EBL in Figures || and 
[To| For instance, if the local UV luminosity density is a 
factor of 5 higher than the value we have adopted and if 
5(2000 A) ~ 1.7 over the range 0 < z < 1, then the rest- 
frame UV luminosity density at z = 1 is similar to that 
measured by CFRS, and the predicted U 300 EBL will be 
roughly 3.5x10“® cgs, very similar to the EBL we derive 



Fig. 12.— K-corrections for various filters as a function of redshift 
calculated using Poggianti (1997) SEDs for present day E, Sa, and 
Sc galaxies. Bandpasses are standard Johnson/Cousins filters and 
the corresponding WFPC2 filters, shown with solid and dotted lines, 
respectively. 


from our modeled local luminosity density and the mean 
values for 5(A) from CFRS. 

4.3. Discussion 

Evolution in the luminosity density of the form (1 -|- 
z)"^^^^ at 0 < z < 1 and slower growth or stabiliza¬ 
tion at z > 1 , such as suggested by redshift surveys at 
0 < z < 4, is consistent with the detected EBL for values 
of 5(A) consistent with CFRS. The strongest constraints 
we can place on the EBL span a factor of 5 in flux. As 
such, stronger evolution between 0 < z < 1 than re¬ 
ported by CFRS or continuing evolution at z > 1 can¬ 
not be tightly constrained. We note that recent results 
from Wright (2000), which constrain the 1.25^m EBL flux 
to be 2 . 1 (±l.l)xl 0 “® cgs, are in good agreement with our 
results, but do not improve the constraints on the high red¬ 
shift optical luminosity density over those discussed above. 

In contrast with our results, previous authors have 
claimed good agreement between the flux in the raw num¬ 
ber counts from the HDF and integrated flux in the mea¬ 
sured CFRS luminosity density to z = 1 under the assump¬ 
tion that the luminosity density drops rapidly at z > 1 
(Madau, Pozzetti, & Dickinson 1998, Pozzetti etal. 1998). 
In that work, the errors in faint galaxy photometry which 
cause ~ 50% underesti mate s of the total light from U > 23 
galaxies (discussed in §4.1) are compensated by the as¬ 


sumption that the luminosity density drops rapidly beyond 
z = 1. That assumption was based on measurements of the 
flux from Lyman-limit systems in the HDF held by Madau 
et al (1996), which are substantially lower than measure¬ 
ments by Steidel etal. (1999) due to underestimates of the 
volume corrections and to the small-area sampling. We 
hnd that the detected EBL is not consistent with luminos¬ 
ity evolution comparable to the CFRS measured values at 
0 < z < 1 if the luminosity density drops rapidly at z > 1 . 


5. FLUX FROM SOURCES BELOW THE SURFACE 
BRIGHTNESS DETECTION LIMIT 
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The fractional EBL23 flux which comes from detected 
sources is simply the ratio of the flux in detected sources 
(measured by “ensemble photometry”) to the detected 
EBL23 (±2cr limits). The maximum fractional EBL23 flux 
coming from undetected sources is what remains: 0-65% at 
C^ 300 ) 0-80% at V 555 , and 0-80% at Isu- Although these 
limits include the possibility of no additional contribution 
from undetected sources, it is worthwhile to note that if 
the progenitor of a normal disk galaxy at z = 0.1 (central 
surface brightness /io(E) ~ 21.3 mag arcsec“^ and E ~ 22 
mag) existed at z ~ 2 , then the z ~ 2 progenitor would 
have a core surface brightness (within 0.2 x 0.2 arcsec^) 
of ^ 26 mag arcsec“^ for standard K-corrections and pas¬ 
sively evolving stellar populations (e.g. Poggianti 1997), 
which is roughly the detection limit for the HDF. In partic¬ 
ular, regardless of the exact evolutionary or k-corrections, 
dimming due to cosmological effects alone (redshift and 
angular resolution) produce ^ 5 mag of surface bright¬ 
ness dimming. This effect is independent of wavelength, 
so that dimming at other bandpasses is similar, modulo 
differences in the evolutionary and K-corrections (shown 
in Figures |^, |l^ and E)- At /, for example, the drop in 
surface brightness for a disk galaxy at z ~ 2 relative to 
z ~ 0 is ~ 0.5 mag greater than for V. The progenitor of 
a typical disk galaxy at z = 0.1, which has E — / ~ 0.9 
and Ho{T) ^ 20.3 mag arcsec“^ (de Jong & Lacey 2000), 
will have E — / ~ 0.5 at z ~ 2 and ~ 25.5. Thus the 
typical disk galaxy at z ^ 2 is close to the HDF detection 
limit in / as well as E. Irrespective of the color evolution 
with redshift, the point is that cosmological surface bright¬ 
ness dimming alone suggests that a significant fraction of 
the EBL23 may come from normal galaxies at redshifts 
z < 4 which are undetectable in the HDF. Furthermore, 
recent redshift surveys for low surface brightness (LSB) 
galaxies now suggest that the distribution of galaxies in 
fio is nearly flat for /ig > 22.05 mag arcsec“^ at some lu¬ 
minosities (Sprayberry etai 1997, Dalcanton etal. 1997, 
O’Neil & Bothun 2000, Blanton etal. 2001, Cross etal. 
2001). If such populations exist at high redshift, they may 
contribute significant flux to the EBL as presently unde¬ 
tectable sources. 

In this section, we explore the possible contributions to 
the EBL23 from galaxies at all redshifts which escape de¬ 
tection in the HDF because of low apparent surface bright¬ 
ness. To do so, we have simulated galaxy populations at 
redshifts 0 < z < 10 as the passively evolving counterparts 
of local galaxy populations and then “observed” the sim¬ 
ulated galaxies through the Gaussian 0.1" FWHM point 
spread function of WFPC2. We define the surface bright¬ 
ness detection threshold to be consistent with the 5cr detec¬ 
tion limits of the HDF images (see Table ||), which corre¬ 
spond to roughly the turn-over magnitude in the number 
counts. This exercise is not meant to approximate real¬ 
istic galaxy populations at high redshift; the evolution of 
galaxy populations in surface brightness, luminosity, and 
number density is so poorly constrained at present that 
more specific modeling is unwarranted. The models dis¬ 
cussed here simply address the question: how much of the 
total flux from a local-type galaxy population at a redshift 
z can be resolved into individual sources? 

5.1. Models 



2468 2468 2468 

redshift 


Fig. 13. — Evolutionary corrections for various filters as a function 
of redshift calculated using Poggianti (1997) SEDs for E, Sa, and Sc 
galaxies evolving passively with redshift. Bandpasses are standard 
Johnson/Cousins filters and the corresponding WFPC2 filters, shown 
with solid and dotted lines, respectively. 
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Fig. 14.— The sum of the evolutionary and K—corrections for John- 
son/Couains (solid lines) and WFPC2 filters (dotted lines) shown in 
Figures and [l3| 


We have considered three models for the central surface 
brightness distributions of disk galaxies in order to explore 
the possible contributions from LSB galaxies. These mod¬ 
els are taken from Ferguson & McGaugh 1995 (FM95) and 
can be generally described as follows: a standard passive 
evolution model in which all galaxies have central sur¬ 
face brightnesses in the range 21 < /io(5j) < 22 mag 
arcsec“^ (Model PE); an LSB-rich model (Model A), in 
which galaxies of all luminosities have 21.5 < fio{Bj) < 25 
mag arcsec“^; and a more conservative LSB model (Model 
B), in which fig is monotonically decreasing for galax¬ 
ies fainter than LjL* < 1 and 21 < p,o{Bj) < 22 mag 
arcsec“^ otherwise (see Figure We include passive lu¬ 
minosity evolution in all three models, as the no-evolution 
models have been clearly ruled out by both number counts 
and our own EBL results. As Model A (Model PE) has a 
broader (narrower) surface brightness distribution than is 
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log (L/L-) 


Fig. 15. — The dark hatched regions in the three panels show the 
surface brightness distribution as a function of luminosity for the 
models adopted here (Model A, Model B, and PE). For comparison, 
the light hatched region shows the surface brightness distribution as 
a function of luminosity (relative to L*) as found by Blanton etal. 
( 2001 ). 



Fig. 16.— For the Johnson/Cousins bandpasses indicated in each 
panel, we plot the EBL from galaxies which are not individually 
detected in each redshift bin, normalized by the total EBL in each 
redshift bin. Models A, B and PE are marked with solid, dashed, 
and dotted lines, respectively. Thick lines correspond to simulations 
run with h = 0.7, flo = 1.0, while thin lines correspond to h = 0.7, 
Qq = 0.2. The dotted vertical line in each panel indicates the Lyman 
limit for the band-pass. 


found by recent LSB surveys (Sprayberry et al. 1997, Dal- 
canton etal. 1997, O’Neil & Bothun 2000), these models 
are taken as illustrative limits on the fraction of low surface 
brightness galaxies in the local universe. Recent determi¬ 
nations of the number density of galaxies as a function of 
both luminosity and surface brightness (c.f. Blanton et al. 
2001 and Cross etal. 2001) are well bracketed by these 
models: Model A allows for too many low surface bright¬ 
ness galaxies, while the PE model clearly allows for too 


few (see Figure |l^ . 

As described in Table 2 of FM95, each surface bright¬ 
ness distribution model has been paired with a tuned lu¬ 
minosity function, so that each model matches the ob¬ 
served morphological distributions and luminosity func¬ 
tions recovered by local redshift surveys. The models 
include identical distributions in the relative number of 
galaxies of different Hubble types (E/SO, SO, Sab, Sbc, 
and Sdm), which are described by bulge-to-disk flux ra¬ 
tios of 1.0, 0.4, 0.3, 0.15, and 0.0, with small scatter. The 
bulge components for all galaxies have E-type SEDs, and 
SO to Sdm galaxies have disk components with E, Sa, Sb, 
and Sc-type SEDs, for which we have used the Poggianti 
(1997) models. Bulges were given r^/’*^-law light profiles 
with central surface brightnesses drawn from the empir¬ 
ical relationship found by Sandage & Perelmuter (1990), 
IJo — —0.48Mbj,-|- 11.02. For disk components, we adopted 
exponential light profiles, with surface brightnesses drawn 
from the 3 model distributions for disk galaxies listed 
above. 

We have calculated passive evolution and K-corrections 
from the population synthesis models and SEDs of Pog¬ 
gianti (1997), shown in Figures |l^ and |^, and we have 
assumed uniform comoving density as a function of red¬ 
shift in all cases. All models were run with Hq = 50 and 
70 km s“^ Mpc“^ and (Om = 0.2, SIa = 0), (1,0), and 
(0.2, 0.8). Different values of Hq have little effect (< 10%) 
on the integrated counts or background. The total back¬ 
ground increases for models with larger volume — (IjO), 
(0.2,0), and (0.2,0.8), in order of increasing volume — but 
the fractional flux as a function of redshift changes by less 
than 10% with cosmological model. 

All three models under-predict the number counts and 
integrated flux in observed sources, as expected, and will 
clearly under-predict the total EBL as illustrated in the 
passive evolution model discussed in §4.2. 


5.2. Results 

In Figures ^ through we plot the distribution of 
the total flux from the modeled galaxy populations as a 
function of redshift, wavelength, and origin (detectable or 
undetectable galaxies). Detection limits applied at each 
bandpass are the 5cr detection limits of the HDF cata¬ 
log (Williams etal. 1996), with appropriate conversions to 
the ground-based filter bandpasses, summarized in Table 
1^ The conversions given in this Table include differences 
in the evolutionary corrections and K-corrections between 
WFPC2 and UBVRI filters (see Figures [T2|-p^, which are 
generally less than 0.3 mag and change by less than 0.1 
mag at z > 0.5. We only consider sources with E > 23 
mag here, and we assume perfect photometry for sources 
which meet the detection criteria. 

In Figure |^, we show the fraction of the total flux 
which comes from undetected sources as a function of red¬ 
shift. For all models, this plot demonstrates that if galaxy 
populations at higher redshifts are the passively evolv¬ 
ing counterparts of those in the local universe, the flux 
from undetected sources becomes significant by redshifts 
of 1 < z < 3. The undetected fraction is the highest in the 
U band, due to the high sky noise and low sensitivity of 
the F300W HDF images relative to the other bandpasses 
which define our detection criteria. The detection fractions 
are similar in B and V, where detection limits and galaxy 
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Table 3 

Adopted Detection Limits of the HDF 


HST 

Filter 

pcore 

AB mag arcsec“2 

m 

AB mag 

Johnson / Cousins 
Filter 

pcore 

mag arcsec~2 

m 

mag 

WFPC2/F300W 

25.0 

27.5 

u 

23.9 

26.4 

WFPC2/F450W 

25.8 

28.3 

B 

25.8 

28.3 

WFPC2/F606W 

26.3 

28.8 

V 

25.9 

28.4 

WFPC2/F606W 

26.3 

28.8 

R 

25.5 

28.0 

WFPC2/F814W 

25.8 

28.3 

I 

25.2 

27.7 



1 3 3 4 5 6 

redshifl 



wavelength {/j.) 


Fig. 17 .— For the Johnson/Cousins bandpasses indicated in each 
panel and the models discussed in the text, we plot the redshift 
distribution of the EBL — the differential EBL from all galaxies as a 
function of redshift, normalized by the total EBL in each band. Line 

g es correspond to the models as described in the caption of Figure 
In this plot, cosmological models are virtually indistinguishable 
ause the fractional volume per redshift bin changes very little 
with n. The dashed vertical line in each plot indicates the redshift 
corresponding to the Lyman limit for the central wavelength of each 
bandpass. 


Fig. 18 .— The lines show the fraction of the EBL that comes 
from undetected galaxies as predicted by our models. Line types are 
as in Figure hrf. Arrows show the upper limits on the fraction of 
the EBL whicn might come froni_ undetected ealaxies based on the 
EBL detections summarized in ^ and Table Q These arrows show 
the ratio of flux recovered by ensemble photometry (from resolved 
galaxies) and the two sigma upper limits of our EBL detection. See 
^.3 for discussion. 


colors are similar. The fraction of light from undetected 
sources in / is small a.t z < 2 due to the generally red 
color of galaxies, but increases beyond that redshift due 
to cosmological effects. Model A, with the largest fraction 
of low /To galaxies, has the sharpest increase in the unde¬ 
tected EBL with redshift, as expected. A balance between 
evolutionary- and K-corrections at 1 < z < 3 slow this 
trend and cause the dip in the fraction of undetected light 
in B, E, and R. The Lyman limit for each band obvi¬ 
ously represents the highest redshift from which one could 
expect to detect flux. 

In Figure |l^, we plot the distribution of light with red¬ 
shift in these models. All three models have roughly the 
same distribution of /ebl(A, z) simply because all models 
employ a uniform comoving number density with redshift 
and the same passive luminosity evolution. Although we 
do not intend to realistically predict the redshift distribu¬ 
tion of the EBL, we show this plot for comparison with 


Figure |l^ to indicate the redshifts at which the majority 
of undetect ed g alaxies lie in these models. Looking at Fig¬ 


ures 


laetected ga 
111 and [Td t' 


_[together, it is clear that while 40— 100% of 

the i?-band flux from z > 3 is in undetectable sources for 
all of the models considered, only a small fraction of the 
total B-band EBL comes from those redshifts. Thus, the 
majority of the light from unresolved sources comes from 
1 < z < 3 at i? for local-type galaxy populations in this 
scenario. 

Figure [T^ shows the fraction of EBL23 coming from un¬ 
detected sources as a function of wavelength. These mod¬ 
els indicate that 10-35% of the light from the high redshift 
counterparts of local galaxy populations would come from 
(individually) undetected sources in bandpasses between 
V and I with sensitivity limits similar to the HDF, 15- 
40% would come from undetected sources at B, and 20- 
70% would come from undetected sources at U. This trend 
with wavelength (smallest fraction of undetected sources 
around 5000A) follows the trend in the dete ction limits 
of the HDF bandpasses, as discussed in §4.1, Note that 
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the color of the EBL23 is similar to the color of detected 
galaxies (see Figure ||) in V and I, as is the 2a lower limit 
of minEBL23 (see also Table ||). 

We stress again that cosmological surface brightness 
dimming and the fraction of LSBs in each model are 
the dominant effects which govern how much light comes 
from undetected sources and these effects are independent 
of wavelength. The passive luminosity evolution correc¬ 
tions, K-corrections, and the HDE-specific detection lim¬ 
its we adopt will determine how the fraction of undetected 
sources varies with wavelength. Finally, we note that al¬ 
though the surface brightness distribution of galaxies as a 
function of redshift is presently unconstrained, and may or 
may not show significant variation with redshift, it is un¬ 
likely that the surface brightness distribution at any red¬ 
shift is significantly more extreme than the distribution 
bracketed by our models. Bearing these uncertainties in 
mind, we can use the results of these models to estimate 
the value of EBL23 based on the minEBL23 (the flux in in¬ 
dividually detected galaxies from ensemble aperture pho¬ 
tometry) and the undetected fractions summarized above. 
If the universe is populated by galaxies with surface bright¬ 
ness distributions like those in the local universe, then 
these models suggest the following values for EBL23: 2.6- 
7.0x10“® cgs, 1.0-1.3x10“® cgs, and 0.9-1.2x10“® cgs at 
U 3 Q 0 , V 555 and / 814 , respectively. These ranges are in good 
agreement with our detected values for EBL23 (see Table 
H), and with the estimates of the EBL2 3 ba sed on the 
corrected number counts we presented in §4.1. 

6. THE BOLOMETRIC EBL (0.1-1000/im) 

In Figure |^, we plot EBL detections to date, together 
with the integrated light in detectable sources (lower lim¬ 
its to the EBL) in units of i/Ii, between 0.1 and lOOO/im.Q 
The DIRBE and FIRAS detections at A > lOO/im and 
the lower limit from IRAS detected galaxies at 10-100/im 
indicate that energy is contained in the far infrared por¬ 
tion of the spectrum. Given that light from stellar nu- 
cleosynthesis is emitted at wavelengths 0.1 — lO/rm, Fig¬ 
ure M emphasizes the fact that 30% or more of the light 
from stellar nucleosynthesis has been redistributed into the 
wavelength range 10 — 1000 /rm by dust absorption and re¬ 
radiation and, to a lesser degree, by cosmological redshift- 
ing. Realistic estimates of the total energy from stellar 
nucleosynthesis must therefore be based on the bolometric 
EBL from the UV to IR. In lieu of accurate measurements 
in the mid-IR range, realistic models of dust obscuration 
and the dust re-emission spectrum (dust temperature) are 
needed. To discuss the optical EBL in the context of star 
formation, we must therefore first estimate the bolometric 
EBL based on the optical EBL detections presented here 
and current measurements in the far-IR. We do so in the 
following section. 

6.1. Models 

Efforts to predict the intensity and spectrum of the EBL 
by Partridge & Peebles (1967) and Harwit (1970) began 

^ The total energy per unit increment of wavelength is given by 
I = J = f ^I\dlnX. By plotting energy as \Ix = vl^, against 
log A, the total energy contained in the spectrum as a function of 
wavelength is proportional to the area under the curve. We give 
\Ix in standard kms units of nW m~^sr“^, equivalent to 10“® ergs 
s ^cm 



Wavelenglh(/im) 

Fig. 19.— EBL detections, limits, and models as a function of 
wavelength. The filed circles show the EBL detections with 2cr error 
bars and lower limit symbols as defined in Figure^ Also plotted are 
lower limits from Armand etal. (1994) at 2000A, the HDF (Williams 
etal. 1996), Gardner etal. (1997) at 2.2/im, IRAS (Hacking &; Soifer 
1991) at lO-100/.im, and Blain et al. (1999) at 450 and 850 fim. These 
lower limits are based on the integrated flux in detected sources at 
each wavelength. Upper limits marked in bold are from Hurwitz, 
Bowyer, &; Martin (1991) at 1600A and DIRBE (Hauser etal. 1998). 
The open circles indicate DIRBE detections. The bold line at 125— 
1000//m shows the FIRAS detection (Fixsen etal. 1998). Detec¬ 
tions in the near-IR are from Wright (2000, filled triangle), Gorjian, 
Wright, Chary (2000, open squares), and Wright & Reese (2000, 
filled square). The lines which indicate models are all labeled and are 
from Malkan Stecker (1998, MS98), Dwek etal. (1^8, D98), and 
Fall, Chariot, Sz Pei (1996, FCP96) as described in The shaded 
region shows the D98 model rescaled to match the range allowed by 
our EBL detections. 


with the intent of constraining cosmology and galaxy evo- 
lution. Tinsley (1977, 1978) developed the first detailed 
models of the EBL explicitly incorporating stellar initial 
mass functions (IMFs), star formation efficiencies, and 
stellar evolution. Most subsequent models of the EBL 
have focused on integrated galaxy luminosity functions 
with redshift-dependent parameterization, with particular 
attention paid to dwarf and low surface brightness galaxies 
(see discussions in Guiderdoni & Rocca-Volmerange 1990, 
Yoshii & Takahara 1988, and Vaisanen 1996). 

More recent efforts have focused on painting a detailed 
picture of star formation history and chemical enrichment 
based on the evolution of resolved sources. The evolu¬ 
tion of the UV luminosity density can be measured di¬ 
rectly from galaxy redshift surveys (e.g., Lilly etal. 1996, 
Treyer etal. 1998, Cowie etal. 1999, Steidel etal. 1999, 
Sullivan etal. 2000), from which the star formation rate 
with redshift, p*{z), can be inferred for an assumed stel¬ 
lar IMF. The mean properties of QSO absorption systems 
with redshift can also be used to infer p*{z), either based 
on the decrease in Hi column density with decreasing red¬ 
shift (under the assumption that the disappearing Hi is 
being converted into stars) or based on the evolution in 
metal abundance for an assumed IMF and corresponding 
metal yield (e.g., Pettini etal. 1994; Lanzetta, Wolfe, & 
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Turnsheck 1995; Pei & Fall 1995). Estimates of the star 
formation rate at high redshift have also come from esti¬ 
mates of the flux required to produce the proximity effect 
around quasars (e.g., Gunn & Peterson 1965, Tinsley 1972, 
Miralda-Escude & Ostriker 1990). Using these constraints, 
the full spectrum of the EBL can then be predicted from 
the integrated flux of the stellar populations over time. 

Unfortunately, all methods for estimating p^,{z) contain 
significant uncertainties. The star formation rate deduced 
from the rest-frame UV luminosity density is very sensitive 
to the fraction of high mass stars in the stellar initial mass 
function (IMF) and can vary by factors of 2-3 depending 
on the value chosen for the low mass cut off (see Leitherer 
1999, Meader 1992). Aside from the large uncertainties in 
the measured UV luminosity density due to incomplete¬ 
ness, resolved sources at high redshift are biased towards 
objects with dense star formation and may therefore paint 
an incomplete picture of the high-x universe. Also, large 
corrections for extinction due to dust must be applied to 
convert an observed UV luminosity density into a star for¬ 
mation rate (Calzetti 1997). 

The SFR inferred from QSO absorption systems, whether 
from consumption of Hi or increasing metal abundance, is 
also subject to a number of uncertainties. In all cases, 
samples may be biased against the systems with the most 
star formation, dust, and metals: dusty foreground ab¬ 
sorbers will obscure background QSOs, making the fore¬ 
ground systems more difficult to study. In addition, large 
scale outflows, a common feature of low-redshift starburst 
galaxies, have recently been identified in the high-redshift 
rapidly star-forming Lyman break galaxies (Pettini et al. 
2000 ), suggesting that changes in the apparent gas and 
metal content of such systems with redshift may not have 
a simple relationship to p*{z) and the metal production 
rate. The mass loss rate in one such galaxy appears to 
be as large as the star formation rate, and the recent ev¬ 
idence for Civ in Ly-a forest systems with very low Hi 
column densities (^xlO^"' cm“^) suggests that dilution of 
metals over large volumes may cause underestimates in the 
apparent star formation rate derived from absorption line 
studies (see Ellison et al. 1999, Pagel 1999, Pettini 1999 
and references therein). 

Finally, regarding the predicted spectrum of the EBL, 
the efforts of Fall, Chariot, & Pei (1996) and Pei, Fall, 
& Hauser (1998) emphasize the need for a realistic distri¬ 
bution of dust temperatures in order to obtain a realistic 
near-IR spectrum. 

With these considerations in mind, we have adopted an 
empirically motivated model of the spectral shape of the 
EBL from Dwek etal. (1998, D98). This model is based 
on p 4 ,{z) as deduced from UV-optical redshift surveys and 
includes explicit corrections for dust extinction and re¬ 
radiation based on empirical estimates of extinction and 
dust temperature distributions at z = 0. The comoving 
luminosity density can then be expressed explicitly as the 
sum of the unattenuated stellar emission, es(;/, z), and the 
dust emission per unit comoving volume, ed(v,z). Equa¬ 
tion 1^ then becomes 



[esjiy, z) + ed(zz, z)\dVc{z) 
A-kDl^zY 


(3) 


D98 estimate the ratio ed(A, 0)/es(A, 0) by comparing the 
UV-optical luminosity functions of optically detected galax¬ 


ies with IR luminosity function of IRAS selected sources. 
Using values of £ = (1.30 ± 0.7) x lO^L© Mpc“^ for 
the local stellar luminosity density at O.l-lO^m and £ = 
0.53 X 1 O®£0 Mpc“^ for the integrated luminosity density 
of IRAS sources, Dwek etal. obtain ed(A, 0)/es(A, 0) ^ 0.3. 
The redshift independent dust opacity is assumed to be 
an average Galactic interstellar extinction law normalized 
at the U-band to match this observed extinction. D98 
then calculate the EBL spectrum using the UV-optical 
observed p*{z), a Salpeter IMF (0.1 < M < 12QMq), stel¬ 
lar evolutionary tracks from Bressan etal. (1993), Kurucz 
stellar atmosphere models for solar metallicity, redshift- 
independent dust extinction, and dust re-emission match¬ 
ing the SED of IRAS galaxies. 

The starting-point UV-optical p*{z) for this model is 
taken from Madau, Pozzetti, & Dickinson (1998), which 
under-predi cts t he detected optical EBL presented in Pa¬ 
per I (see §0). While D98 discuss two models which 
include additional star formation at z ^ 1, the additional 
mass is all in the form of massive stars which radiate in¬ 
stantaneously and are entirely dust-obscured, resulting in 
an ad hoc boost to the far IR-EBL. We instead simply 
scale the initial Dwek etal. model by x2.2 to match the 
2cr lower limit of our EBL detections and x4.7 to match 
the 2 a upper limit, in order to preserve the consistency of 
the D98 model with the observed spectral energy density 
at z = 0. In that any emission from z > 1 will have a red¬ 
der spectrum than the mean EBL, simply scaling in this 
way will produce a s pectrum which is too blue. However, 
as discussed in §0, it is also possible that the z < 0.5 
UV luminosity density has been underestimated by opti¬ 
cal surveys, so that the bluer spectrum we have adopted 
may be appropriate. Note that the resulting model is in 
excellent agreement with recent near-IR results at 2.2 and 
3.5/im (Wright & Reese 2000; Gorjian, Wright, & Chary 
2000; Wright 2000) and also with the DIRBE and FIRAS 
results in the far-IR. Adopting this model, we estimate 


where errors are Icr errors associated with the fit of that 
template to the data. 

Due to the corrections which account for the redistribu¬ 
tion by dust of energy into the IR portion of the EBL, the 
star formation rate implied by the unsealed (or scaled) 
D98 model is 1.5 (or 3.3-7.1) times larger than the star 
formation rate adopted by Madau etal. (1998). The dust 
corrections used by Steidel et al. (1999) produce a star for¬ 
mation rate which is roughly 3 times larger than used in 
the unsealed D98 model, slightly smaller than the scal¬ 
ing range adopted here, which is consistent with the fact 
that the CFRS and Steidel etal. (1999) luminosity densi¬ 
ties are slightly bel ow our minimum values for the EBL, 


as discussed in §4.2 


6.2. Energy from Accretion 

As mentioned briefly in §|^, another significant source 
of energy at UV to far-IR wavelengths is accretion onto 
black holes in AGN and quasars. The total bolometric 
flux from accretion can be estimated from the local mass 
function of black holes at the centers of galaxies for an 
assumed radiation efficiency and total accreted mass. Re¬ 
cent surveys find Mbh ~ 0.005Msph, in which Mgph is the 
mass of the surrounding spheroid and Mbh is the mass of 
the central black hole (Richstone etal. 1998, Magorrian 
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etal. 1998, Salucci etal. 1999, and van der Marel 1999). 
Following Fabian & Iwasawa (1999), the energy density in 
the universe from accretion is given by 

fbh = 0.0050,phPcritc2 7-^7^ , (4) 

(1 + ^e) 

in which T^bh is the radiation efficiency, flsph is the observed 
mass density in spheroids in units of the critical density, 
Pcrit, and (1 + Ze) compensates for the energy lost due 
to cosmic expansion since the emission redshift Zg. The 
bolometric flux from accretion is then 


z = 0 in Equation ^ then simplifies to 

Thai _ ^ 

47 r(l + Ze) 


( 6 ) 


in which £* is the bolometric energy density from stellar 
nucleosynthesis and (l + Zg) compensates for energy lost to 
cosmic expansion. In the case of instantaneous formation 
and emission, £* can be expressed in terms of the total 
energy released in the nucleosynthesis of He and heavier 
elements: 

£* = £l*PcritC^v{^Y + Z) , (7) 


f 77^ 10 ^ nWm-^sr-i (5) 

47r (1 + Ze ) 

for rybh ~ 0.1, Ze ~ 2, pcrit = 2.775 x 10“ M© Mpc-^, 
Ho = lOO/i km s“^ Mpc“^, and Hsph ~ 0.0018to°ooo85^~^ 
(Fukugita, Hogan, & Peebles 1998, FHP98). 

The observed X-ray background (0.1-60 keV) is ~ 0.2 nW 
m“^ s“^. The large discrepancy between the detected X- 
ray flux and the estimated flux from accretion has led to 
suggestions that 85% of the energy estimated to be gen¬ 
erated from accretion takes place in dust-obscured AGN 
and is emitted in the thermal IR (see discussions in Fabian 

1999) . Further support for this view comes from the fact 
that most of the soft X-ray background (below 2 keV) 
is resolved into unobscured sources (i.e., optically bright 
quasars), while most of the hard X-ray background is as¬ 
sociated with highly obscured sources (Mushotzky etal. 

2000) . Photoelectric absorption can naturally account for 
the selective obscuration of the soft X-ray spectrum. Best 
estimates for the fraction of the far-IR EBL which can 
be attributed to AGN are then < lOh nW m“^ sr“^, or 
< 30% of the observed IR EBL. This is in good agreement 
with estimates of the flux from a growing central black hole 
relative to the flux from stars in the spheroid based on ar¬ 
guments for termination of both black hole accretion and 
star formation through wind-driven ejection of cool gas in 
the spheroid (Silk & Rees 1998, Fabian 1999, Blandford 
1999). Together, these studies suggest that < 15% of the 
bolometric EBL comes from accretion onto central black 
holes. 


7. STELLAR NUCLEOSYNTHESIS: ft, AND \Z\. 


The bolometric flux of the EBL derived in §S.l 


IS a 

record of the total energy produced in stellar nucleosyn¬ 
thesis in the universe, and so can be used to constrain 
estimates of the baryonic mass which has been processed 
through stars. The relationship between processed mass 
and background flux depends strongly on the redshift de¬ 
pendence of star formation and on the stellar IMF, but is 
only weakly depende nt o n the assumed cosmology for the 
reasons discussed in §4.2 and §||. 

As an illustrative case, we can obtain a simple estimate 
of the total mass processed by stars by assuming that all 
stars formed in a single burst at an effective redshift Ze, 
and that all the energy from that burst was emitted in¬ 
stantaneously. The assumption of instantaneous emission 
does not strongly affect the result because most of the 
light from a stellar population is emitted by hot, short¬ 
lived stars in the first ~ lOMyr. The integrated EBL at 


in which 77 (^ 0.0075) is the mean conversion efficiency 
of energy released in nuclear reactions and AY and Z are 
the mass fractions of ^He and metals. Inverting Equation 
1 ^ the total baryonic mass processed through stars in this 
model can be derived from a measurement of the bolomet¬ 
ric EBL using the expression: 


4^(1+ ze) 

^rjpcrit {AY + Z) 


We can bracket a reasonable range for {AY + Z) by 
assuming the solar value as a lower limit, and the mass 
weighted average of the metal conversion fraction in E/SO 
and spirals galaxies as the upper limitAssuming a 3:2 
ratio of E/SO to Sabc galaxies (Persic & Salucci 1992), 
we find {AY + Z) = 0.25 ± 0.15. For Zg = 1.5, the to¬ 
tal baryonic mass processed through stars corresponding 
to a bolometric EBL of 100 ± 20 nW m“^sr“^ is then 
Hh. = 0.0030(±0.0019)/i“^ in units of the critical density, 
or 0.16(±0.10)nB for Hb = 0.019(±0.001)/i-2 (Buries & 
Tytler 1998). Again, this calculation assumes a single red¬ 
shift for star formation with all energy radiated instanta¬ 
neously at the redshift of formation. 

The true history of star formation is obviously quite dif¬ 
ferent from this illustrative case. For time-dependent emis¬ 
sion and formation, the bolometric EBL is the integral of 
the comoving luminosity density corresponding to realis¬ 
tic age- and redshift-dependent emission (see Equation ^. 
For comparison, instantaneous formation at the same red¬ 
shift assumed above {zg = 1.5) with a modified Salpeter 
IMF and time-dependent emission based on SEDs from 
Buzzoni (1995) would imply = 0.0037(±0.0007)ft.“^ for 
our estimate of the bolometric EBL (for details see MPD98 
and Madau & Pozzetti 1999, MP99). The mean of this 
estimate is about 20 % higher than that from the instanta¬ 
neous formation and emission model discussed above. The 
two models are very similar because the vast majority of 
energy from a stellar population is emitted in the first 
~ lOMyrs. The quoted uncertainty is smaller than for our 
illustrative model because the error range reflects only the 
uncertainty in our estimate of the bolometric EBL and no 
uncertainties in the adopted IMF. 

® Solar values of AY and Z are 0.04 and 0.02, implying DYjDZ = 
2. Interstellar absorption measurements of DYjDZ in the solar 
neighborhood are closer to the range 3-4, implying AY ~ 0.07. He¬ 
lium white dwarfs contribute an additional 10% of the local stellar 
mass to the estimate of AY (Fleming, Liebert &; Green 1986), so that 
we have (AY + Z) = 0.07 + 0.02 -|- 0.1 ~ 0.2 as a local estimate for 
systems with solar metallicity. This is similar to estimates for other 
local spiral galaxies. Estimates for E/SO galaxies are as high as 0.5 
(Pagel 1997). (Note that the He mass produced in stars is written 
as AY to distinguish it from the total He mass, which includes a 
primordial component.) 
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For the same IMF and SFDs, a redshift-dependent star 
formation rate for 0 < z < 4 based on the observed UV 
luminosity density and taking dust obscuration into ac¬ 
count (see Steidel etal. 1999) would imply that almost 
twice as much mass is processed through stars than in the 
instantaneous formation model above (MP99). Relative to 
the instantaneous-formation models, the same bolometric 
FBL flux corresponds to a larger value of D* when we con¬ 
sider time-dependent star formation because more of the 
emission occurs at higher redshifts, resulting in greater en¬ 
ergy losses to cosmic expansion. For our estimate of /ebl 
and the calculations of MP99 discussed above, we there¬ 
fore estimate that total mass fraction processed through 
stars is D* = 0.0062(±0.0012)/i-2 or 0.33 ± 0.07Db. We 
adopt this value for the remainder of the paper. 

For this value of the total processed mass, we can cal¬ 
culate the corresponding metal mass which is produced in 
stellar nucleosynthesis. To do so requires an estimate of 
the metal yield — the mass fraction of metals returned to 
the ISM relative to the mass remaining in stars and stellar 
remnants. Best estimates of the metal yield, yz-, lie be¬ 
tween 0.01 (corresponding to a Scalo IMF) and 0.034 (as 
observed in the Galactic bulge) (Pagel 1987, 1999). These 
values incorporate the full range predicted by IMF models 
and observations (see Woosley & Weaver 1995; Tsujimoto 
etal. 1995; Pagel & Tautviaisiene 1997; Pagel 1997). For 
Z 0 = 0.017 (Grevesse, Noels, & Sauval 1996), this metal 
yield range in solar units is = 1.3 ± 0.7^0. If the mass 
fraction remaining in stars and stellar remnants is /, then 
the predicted metal mass density is given by 

= fyz^* , (9) 


and is, obviously, less than the upper limit — the total 
baryon mass fraction from Big Bang nucleosynthesis and 
deuterium measurements. 

We can also compare the metal mass fraction predicted 
by the EBL with the observed mass fraction in gas, stars, 
and stellar remnants in the local universe. Based on re¬ 
cent estimates by FHP99, ^ 80% of the observed baryons 
at z ~ 0 are located in the intracluster gas of groups 
and clusters, 17% are in stars and stellar remnants, and 
only 3% are in neutral atomic and molecular gas. The 
observed metal mass fraction in hot intracluster gas has 
been estimated to be at least 0.33^0 for rich galaxy clus¬ 
ters and 0.25 — \Zq in groups based on X-ray observa¬ 
tions of Fe features (Renzini 1997, Mushotzky & Loewen- 
stein 1997). More recent estimates from Buote (1999, 
2000) based on more detailed models of the tempera¬ 
ture distribution of the intracluster gas have found val¬ 
ues closer to IZq in several clusters and elliptical galax¬ 
ies. For a total cluster and group gas mass density of 
flgas = 0-011-o’oo5^~^ (FHP98) and assuming a metal 
mass fraction of 0.65 ± 0.35^0 for clusters and groups of 
all masses, the observed metal mass fraction in clusters is 
^^Z,gaS = 0-007)1 q QQg/l Zq. Repeating this exercise for 
the stellar component, we assume that the metallicity of 
stars at z ^ 0 is roughly solar (1.0 ± O. 25 Z 0 ) and that 
the mass density in stars is Dgtars = 0.0025(±0.001)/i“^. 
The total metal mass in stars and stellar remnants lo¬ 
cally is then — O.OO25(±O.OOl)/i“^Z0. The to¬ 

tal metal mass fraction in the local universe is then U,z = 
flz.star+Oz,gas = O.OOObtHo^”^ 0Z ■ 
This estimate is consistent with the value indicated by the 
bolometric EBL, O.24(±O.13)Z0 Db, as calculated above. 


which gives ilz = 0.0040(±0.0022)/i ‘^Zq, or flz = 
O.24(±O.13)Z0Db, for an assumed lock-up fraction of 
/ = 0.5. 


Note that we have assumed that the full flux of the EBL 
is due to stellar nucleosynthesis in the above calculations 
of D* and ^lz■ If ^ lO/i n W in ~^sr~^ of the IR EBL is due 
to AGN, as estimated in (36.2, then the energy emitted by 


stars is smaller by about about 7%, and the inferred mass 
fractions are then smaller by about 7% as well. 


7.1. Comparison with other Observations 

The total mass processed by stars is not a directly 
observable quantity because some fraction of the pro¬ 
cessed mass will be hidden in stellar remnants or re¬ 
cycled back into the ISM. Estimates of recycling frac¬ 
tion range from 30-50% for various IMF models (see 
discussions in Pagel 1997), but the cumulative effect of 
many generations of star formation and repeated recy¬ 
cling is difficult to estimate. Firm lower and upper lim¬ 
its for D* are, however, directly observable: the observed 
mass in stars and stellar remnants at z ^ 0 is a lower 
limit to the total mass which has been processed through 
stars, and the total baryon fraction from Big Bang nucle¬ 
osynthesis is an upper limit. FHP98 estimate the mass 
fraction in stars and stellar remnants at z ~ 0 to be 
llstars = 0.0025(±0.001)/i“^, corresponding to a mass- 
to-light ratio of {M/L)b ~ 5.9{M/L)q. In units of Db, 
this lower limit is Dgtars = 0.13(±0.05)hDB. Our esti¬ 
mate of the total mass fraction processed through stars, 
Dh. = 0.33(±0.07)Db, is comfortably above this lower limit 


8. SUMMARY AND CONCLUSIONS 

Based on surface photometry from HST/WFPG2 and si¬ 
multaneous ground-based surface spectrophotometry from 
Las Campanas Observatory, we And mean values for the 
flux of the EBL23 (the background light from sources 
fainter than V = 23ABmag) as follows: Ifsoow = 4.0 ± 
2.5, /f 555 W = 2.7 ± 1.4, and /f 814 W = 2.2 ± 1.0 in 
units of xl0“® cgs, where uncertainties quoted are Itr 
combined statistical and systematic errors. These results 
are presented in detail in Paper I and are summarized 
in §pl Adding in the flux from sources brighter than 
E = 23 ABmag (see Table ^), we find the total EBL 
flux is /f 300 W = 4.3 i 2.6, /fsssw = 3-2 i 1-5, and 
7 f814w = 2.9 ± 1.1 cgs. 

In the context of these measurements of the EBL, we 
have discussed constraints on the slope of number counts, 
the luminosity density as a function of redshift, the frac¬ 
tion of galaxies which lie below current surface brightness 
detection limits, and the history of stellar nucleosynthe¬ 
sis and metal production in the universe. We reach the 
following principle conclusions: 

(1) We find that the corrected number counts at V and 
I magnitudes fainter than 23 AB mag obey the relation 
N oc 10“”" with a = 0.33 ± 0.01, and a = 0.34 ± 0.01, 
respectively, which is consistent with the slope found at 
brighter magnitudes (e.g. Smail etal. 1995, Tyson 1988). 
This is significantly steeper than the slope of the raw HDF 
number counts (~ 0.24 ± 0.1 at E >23 ABmag, and 
^ 0.22± 0.1 at I > 23 ABmag). In contrast with the raw 
counts, the corrected counts show no decrease in slope to 
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the detection limit. If we integrate the corrected num¬ 
ber counts down to an apparent magnitude corresponding 
roughly to a dwarf galaxy (My ~ —10 mag) at z ~ 3, 
y ~ 38 AB mag, we obtain a total flux of 1.2x 10“® cgs in 
both V and I. This is 1.2ct below the mean EBL23 flux 
we estimate at Vgoe (.^ 814 )) suggesting that number counts 
would need to be steeper over some range in apparent mag¬ 
nitude fainter than the current detection limits in order to 
obtain the mean EBL flux we detect, or that the value of 
EBL23 is roughly ~ Icr below our mean detections at V 
and I. 

(2) Based on a local luminosity density consistent with 
Loveday etal. (1992), passive evolution in the luminos¬ 
ity density of galaxies under-predicts the EBL by factors 
of roughly 3, 2, and 2 at Usoo, V 555 , and / 814 , respec¬ 
tively. Note, however, that if the local luminosity density 
is a factor of two higher than the Loveday et al. values we 
have adopted here, as found by Blanton etal. (2001), then 
passive evolution agrees with the flux in resolved galaxies 
(minEBL23) and with our mean EBL detections to within 
1 ( 7 . The mean detected EBL therefore requires stronger 
evolution in the luminosity density than passive evolution 
will produce, however, the exact form of that evolution is 
not well constrained by our results. 

Adopting the local luminosity density assumed by Lilly 
etal. (1996, CFRS), the Icr upper limits of the cumulative 
flux measured by Lilly et al. from redshifts 0 < z < 1 is 
smaller than the flux in resolved sources by more than a 
factor of 2 : this fact alone demonstrates that significant 
flux must be contributed by galaxies at redshifts z > 1. If 
we adopt £(A, z) cx (1 -I- z)'^^'') for the luminosity density 
at 0 < z < 1 based on the Lilly et al. results, then constant 
luminosity density at z > 1, such as suggested by Steidel 
etal. (1999) is consistent with the detected flux in sources 
at V 555 and / 814 , and with the detected EBL at Usoo- At 
the upper limit of the EBL detections, we find that the 
luminosity density can continue to rise as a power law to 
z ~ 2.5 without over-predicting the EBL. 

(3) We have modeled the effects of cosmological K-corrections, 
passive evolution, and (l-|-z)^ cosmological surface bright¬ 
ness dimming on the detectability of local-type galaxy 
populations as a function of redshift. For these mod¬ 
els, we have adopted the spatial resolution and surface 
brightness limits of the HDF. For models which bracket 
the observed surface brightness distribution of galaxies in 
the local universe, we find that roughly 10-40% of the 
EBL from galaxies fainter than E ~ 23 (i.e. those sam¬ 
pled in an HDF-sized image), comes from galaxies which 
are, at present, individually undetectable at wavelengths 

A > 45 OOA, and roughly 20-70% comes from individually 
undetected galaxies at A < 4500A. Most of the flux from a 
local-type galaxy population located at z = 3 would come 
from sources that would not be individually detected in 
the HDF. Our models indicate that the true EBL is likely 
to be between the mean detected EBL23 values and the 
1(7 lower limits of those detections at V and /, and within 
±1(7 at U. 

(4) Scaling the model of the bolometric EBL derived by 
Dwek etal. (1998), which is based on a combined UV- 
optical estimate of the star formation rate and a model for 
dust obscuration and re-emission based on the spectrum 
of IRAS sources, we find that the optical EBL we detect 
corresponds to a total bolometric EBL (0.1 to lOOO/rm) of 


100±20 nW m-^sr-i. 

(5) From this estimate of the total bolometric EBL, we 
estimate that the total baryonic mass processed through 
stars is = 0.0062(±0.0012)/i-2 = 0.33(±0.07)Hb, and 
that the mean metal mass density in the universe is 
Hz = O.OO4O(±O.OO22 )/i-2Z0 = O.24(±O.13)%0Hb, for 
Db = 0.019(±0.001)/i-2 (Buries & Tytler 1998). These 
estimates are consistent with limits from other observa¬ 
tional constraints. 
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